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We investigate the utility of parity detection to achieve Heisenberg-iimited phase estimation for 
optical interferometry. We consider the parity detection with several input states that have been 
shown to exhibit sub shot-noise interferometry with their respective detection schemes. We show 
that with parity detection, all these states achieve the sub-shot noise limited phase estimate. Thus 
making the parity detection a unified detection strategy for quantum optical metrology. We also 
consider quantum states that are a combination of a NOON states and a dual-Fock state, which 
gives a great deal of freedom in the preparation of the input state, and is found to surpass the 
shot-noise limit. 
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I. INTRODUCTION 

Quantum optical metrology deals with the estimation 
of an unknown phase by exploiting the quantum na- 
ture of the input state under consideration Due 
to the inherent uncertainty imposed by quantum me- 
chanics, the problem reduces to minimizing the uncer- 
tainty of the expectation value of a suitable observ- 
able laaiiEiiTiiHiiiiiniiiiiiiiiiiiiTiiisiiTniiiTi 

[181 [191 EHl Ell Ea ESI Ei [23 ES]. in the usual classical 
setting, for a given mean number of input photons, N , the 
phase estimate scales as I/a/ZV, which is usually referred 
as shot-noise limit (SL). It has been shown that by ex- 
ploiting the signature quantum properties such as entan- 
glement the uncertainty can be reduced to the Heisenberg 
limit (HL) of 1/-/V; an enhancement of a factor of \/iV. 
Achieving a sub-shot noise limit or the Heisenberg limit 
depends on the nature of the input states and the detec- 
tion strategy of the output measurement [T31 [THl 1^ . 

Precise optical phase measurement has been an open 
problem for many years and has many applications, most 
notably the gravitational wave detection [27l [28] . Vari- 
ous input states of light and measurement schemes have 
been shown to surpass the SL. Sanders and Milburn have 
computed the optimal measurement |19j . written as Pos- 
itive Operator Valued Measure (POVM), to achieve the 
HL. The method specified by Sanders and Milburn is in- 
dependent of the system phase and thus an optimal one. 
Berry and Wiseman have considered the optimal POVM 
and derived the optimal input state to achieve the min- 
imum uncertainty that scales as the HL [22 . They also 
showed a way of approximately implementing the optimal 
POVM using a feedback technique. With this feedback 
technique along with the Kitaev algorithm for phase esti- 
mation |:42;, Higgins et al. have experimentally achieved 
the HL scaling of the optical phase measurement [l^ . 

In this article, we will discuss the parity measurement, 
which detects whether the number of photons in a given 
output mode is even or odd. The scope of the article 
is restricted to the local phase estimation as opposed 
to the global [TH EU], such that the phase estimation 



can be evaluated by using the linear error propagation 
formula [4 1. The primary objective of this paper is to 
investigate the utility of parity detection scheme and 
not the role of the quantum correlations-also referred 
as entaglement-in achieving the sub-shot -noise limited 
phase estimate. The precise question on the role of en- 
tanglement in achieving sub-shot-noise limited interfer- 
ometry is beyond the scope of this article. The curious 
reader may find Ref. [57] useful in that direction. 

As we make extensive use of the Schwinger repre- 
sentation to analyze the Mach Zehender Interferometer 
(MZI), we wish to present a brief discussion of the rep- 
resentation. Any four-port optical lossless device, such 
as the MZI considered here, can be conveniently de- 
scribed using the Schwinger representation of the an- 
gular momentum |15j . The operators, which form an 
SU(2) rotation group, and describe the MZI jUj are: 
4 = {a^b+ab^)/2jy = {a)b-ab^)/2ij,_ = {a^a-b^b)/2, 
where a and b are the mode operators which obey bosonic 
commutation relation, [d, a^] = [b,b^ — 1. The angular 
momentum operators obey [Ji,Jj] = ietjkJk- The total 
photon number is N ^ a) a+wb, and = J^ + Jy+J^ — 
{N /2){N /2 + 1) is the Casimir invariant. The generator 
for beam-splitter transformation is usually represented 
by J'x [19,1 21 . The combined two mode input state is rep- 
resented by the simultaneous eigen state of and Jz, i.e 
\jlJ-) zi where \j+fj) and {j—fj) represent \N)a and |A^)f, re- 
spectively and j = A^/2 for a fixed input photon number 
N. Correspondingly, if |jV)z represents the output state, 
then 1/ represents the output photon number difference, 
{Na ~ Nii)\out- In this representation the Mach Zehender 
Interferometer is given by an operator, I — exp(—iipjy). 
For a given input state | j/j,) , the output state can be writ- 
ten as, i |j» = e-'^-'y = EL-j IJ^), where 
dl ^{(p) is the usual rotation matrix elements [Ml H5] . 

The paper is organized as follows. In section|TT]we dis- 
cuss the parity measurement and setup a general frame- 
work of calculating the expectation value for an arbitrary 
input state. In section [Hi] we apply it to specific input 
states including a combination of a NOON state and a 
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dual-Fock state, followed by section IV with conclusions 




FIG. 1: Schematic of a Mach-Zehender Interferometer (MZI). 
\tp)in represents the joint input state at ai„ and bi„.The pho- 
ton number states entering at ai„, (bin) are |A'')a(|A'')i,). The 
symbol ip represents the relative phase between the modes, 
within the interferometer. BS: Beam Splitter. 



II. PARITY DETECTION 

Parity detection was first proposed by Bollinger et al. 
in 1996 to study spectroscopy with a "maximally entan- 
gled" state of trapped ions [4]. The detection considered 

there is (—1)^=, where \N)^ represents N atoms in ex- 
cited state. It is straightforward to draw parallel between 
the two level atom and the MZI depicted in Fig. [T] In the 
case of the optical phase measurements, we have {ip) n the 
NOON state. Formally, in the present context of phase 
measurement, considering the MZI shown in Fig. [T| the 
NOON states are written asEilSnilSlIMlEllSSlIMlET]: 



IN 



\N)a>\0)b' + \0)a'\N),, 
V2 



(1) 



where N is the number of photons and a' and b' are the 
two internal modes of the MZI. 

The detection operator which was first proposed by 
Gerry 



P =(-!)' 



(-1) 



3 -J. 



(2) 



There is no particular reason to perform such a detec- 
tion on bout and choosing Gout works equally well. Gerry 
and Campos have applied this operation to interferome- 
try with the NOON states, resulting in the exact HL [?T] . 
As stated in the introduction, these states are not the in- 
put states of the MZI, but are states after the first beam- 
splitter (BS). One can quite easily write down the in- 
put state using a beam-splitter transformation as we will 



do it now. In the Schwinger representation, the NOON 
states is represented as: = + \j,—j))/V2. 

Then we get the actual input state: 



|V.), = e-^t-^^|V')^= E ^mIJ-m), (3) 

where the coefficient is given by 
1 



We used e * 2 



4 ; Jx — gj 5 g I ; g 1 5 jji obtaining the 



above result. Thus 1?/')^ is the input state of the MZI to 
get the exact HL with parity detection. The real part of 
the coefficients are plotted in Fig. [2] 




FIG. 2: The coefficients of the state given in Eq. ^ for an 
input photon number, 7V=100. This is the input state to the 
MZI shown in Fig.[T]such that it is NOON state after the first 
beam splitter and thereby giving the Heisenberg-limited phase 
estimate. 

Recently, Uys and Meystre [TT] noted that a state 
whose coefficients look alike the coefficients plotted in 
Fig. |2] also gives an exact Heisenberg limit. They ob- 
tained this result via numerical simulations and an ex- 
plicit mathematical expression of the input state was not 
given. From Fig. [2] it appears that the state obtained in 
Ref. [TT is a beam-splitter transformation of the NOON 
state which is given by Eq. (|3|. 

For a local phase estimation the phase uncertainty is 
typically given as [4T| : 



6ip 



AP 

K{P)\" 



(4) 



where AP = y (P^) - (P)^ = \/l - {P^, since P^ = 1. 

The usual approach in phase estimation one wishes to 
take, is to fix the detection strategy, i.e. to fix a partic- 
ular phase dependent observable at the output and look 
for the behavior of the input states. While a detection 
of the signal is done at the output ports, one can al- 
ways transform the observable through the second beam 
splitter and can think of the detection within the interfer- 
ometer. Obviously, both descriptions are equivalent but 
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the latter may help understanding the direct effect of 
the measurement on the states after acquiring the phase 
shift. 

We can obtain more insight on P by transforming 
through the BS and further express it as a projection 
operator by using the completeness relation: 

= E (-l)-''<^(-'r)b-^)0-A*|. (5) 



used 



-iipjy 



Ijm) 

The summation 



over 2j describes the situation where the number of 
photons are not fixed. In what follows we will use a 
fixed number of input photons N, and this summation 
will then be dropped. 



III. APPLICATION OF THE PARITY 
DETECTION 



It is straight forward to see that = 1. In terms 
of the photon numbers, and noting that dl ^(—n) = 
we get m 



N 



(6) 



fe=0 



If the NOON states are under consideration for the 
above expression of the observable Q, the only rele- 
vant terms would be for fc = and k = N. Thus 
the observable considered in Refs. [TSJ [SSJ ISHl HO] is: 
Qn = |0,7V) {N,0\ + \N,0) {0,N\, that leads to the HL 
for NOON states. This gives the same expectation value, 
and thus Qn can be implemented with parity detection 
for the NOON states. Explicitly, after the phase shifter, 
\^p)j^ becomes: = |iV,0) + e^^^ |0, A^), which leads 

to (see also Ref. fH]): 



^+isinAfV3, N odd. 



even. 



(7) 



Using Eq. ^ and Eq. Q we immediately get 6ip — 1/N. 

Now we shall obtain the expectation value of parity 
observable for an arbitrary input state. In the Schwinger 
notation an arbitrary two-mode input and the corre- 
sponding output states are written as, 

oo j 

l^)in= E E V-^^^.IJ-M), (8) 
oo j 

Wout = ^ l^>in = E E \J, t^) , (9) 

where ipf^j is the amplitude for the arbitrary input state. 
Using the above equation it is straightforward to calcu- 
late the expectation value of P for an arbitrary input 
state 

(P)out= out(V'l^'IV')out 
oo j 

= E E i-iy-'''r,',,^,.K',(^^)- (lo) 

2i=0 fi,ti' = -j 



A. Parity detection with uncorrelated states 

We begin with a coherent state input at mode ai„ . We 
have at the input: 



2j=0 



7(27)1 



(11) 



where \a\ = n, the average photon number. Using Eq. 
p^, we have: 



2j=0 



exp 



(2j) 



|a|Vl + cos(2y)) 



(12) 



which, in the limit — *■ 0, according to Eq. Q, im- 
mediately leads to 6(pa — l/V^- We thus recover the 
shot- noise limit. This can also be obtained by mea- 
surement at the output, which corresponds to the photon 
number difference. However, in the case of a Jz measure- 
ment, the shot-noise limit is reached when (p tends to odd 
multiples of 7r/2 [IT. 

The next simplest uncorrelated state is a number state 
at Clin and vacuum at bin- Thus the input state: — 

W)a\^)b = J2i=-j^i^,j\jJ)- The subscript s denotes 
the input at a single port of MZI. Using Eq. (10 1, we 
obtain: 



(P)s = 



[l+cos(2y)]^ 
2J 



(13) 



for which, in the limit ip 0, we get dips ~ l/\/2j ~ 
1/^/N. This result shows that parity detection gives the 
same result as the Jz measurement for a single-port input 
state [47] . 

Now we consider the dual-Fock input state [3J M 1211 
I25j . Campos et al. have shown the utility of the par- 
ity measurements for the dual Fock input state |48| . 
Their analysis also includes a comparison and contrast 
of the quantum state distribution with the interferom- 
eter with the NOON states (as these are the states 
within the interferometer). The dual Fock state can be 
In obtaining the above result, we written as: IV')^ = 1-^,-^) = S^=-j '^m,o Ijm)- Using 



4 




FIG. 3: The variance of 5ip versus A'', the input number of photons. (I) The dual fock state, shown as the continuous curve which 
is same as given in Ref. [48]. (II) For the combined state |V')ci in ths absence of relative phase, where (a) a = y^, (i — \J\, d = 0, 

(b) Q = \J\,P= and 6* = 0. (Ill) For a fixed a = /? = ^, (a) £» = 0, (b) represented by x for 61 = tt, (c) 6^ = f . In (III) 
the plots for the case (a) and (b), almost overlap. In all plots the dot-dashed curve is the shot-noise limit (SL) and the dotted 
curve represents the Heisenberg limit (HL). 



Eq. (10) we immediately get {P)d = (-1)^^;^ o(2(^). Us 



ing {P)d i n Eq. ([4|) for a small phase, ^ 0, we get 
5ipd = l/v/2j(i+ 1) cx l/N (see also Ref. [46]). We plot 
this in Fig. [sj^I), which is the same as shown in Ref. [15] . 



B. Parity detection with a combined state 

As discussed above, both the dual- Fock state and the 
NOON state, have an Heisenberg limited phase variance 
with the parity detection. 

It is natural to ask how precisely a state has to be 
prepared to take the advantage of parity or, in general, 
any detection scheme. We now attempt to answer this 
question by considering a combination of a NOON state 
[see Eq. (Is])] and a dual-Fock state such as 



— Cjv 01 



e 2 



V2 



-f3\j,0) 



(14) 



Writing a = |a|e*^° and /3 = |/3|e*^'', where the normal- 
ization constant is given by 

Cat = [I + 2\/2|a||/3|4o(7r/2) cos{0 - (15) 

where 9 — 9^ — Op is the relative phase. With such a 
quantum state as the input of the MZI in Fig. [l] we have 
the output |V')c|out = ^ l'^)c- It gives rise to: 

(V'l ^'IV')c|out = C^{|ap(P)NOON + |/3p(P)dual 

+iJ2%/2|a||/3|d^oCos(7V^)cos(6i)}, (16) 

where 

(P)noon = (-l)^'(e^^'^ + (-I)^e-^^'^)/2, 
(P)duai = (-l)^d^o,o(2^)- 



The above result is obtained using the symmetric proper- 
ties of the rotation matrix elements [44] and the Baker- 
Campbell-Hausdorff (BCH) formula [JT. We can now 
use this in Eq. ^ to calculate Sip. To gain more insight, 
we will first take 9 — Q and look for the Sip for various 
combinations of \a\ and using jap -I- = I. The 
results are shown in the top portion of Fig. Sj^II). 

Next, we fix \a\ = |/3| = 1/^2 and vary 9. This is 
plotted in the in the bottom portion of Fig. IsFlII). In 
both the cases we can clearly see that with at the 
input, using parity detection, we do get the sub-shot noise 
limited variance of the optical phase. This result implies 
a wider applicability of the parity detection. 



C. Parity detection with correlated states 

In this section we consider the quantum-correlated 
states that have been already shown to achieve sub-shot 
noise limited phase estimate under several different de- 
tection schemes. Here we show that all these states does 
achieve the goal with the parity detection and thus mak- 
ing the parity detection, a unified detection scheme. 

Let us begin with the Yurke state. The Yurke state is 
formally written as [\N/2)a\N/2)i, + \{N + l)/2)a\{N - 
l)/2){,]/-\/2. Yurke et al. [T^ have shown that by using 
the output photon difference, a minimum phase sensitiv- 
ity of 2/N can be achieved. In the Schwinger notation 

Yurke state is = E^^-j 73 ((5^,0 + (^^,1) IJ'M)- Using 



Eq.(lO) we get the expectation value of parity as [46] 



{P)y = 



iU2^)-di,{2^) + 2dU2^)]. (17) 



Again in the limit 1^9 ^ 0, we have using Eq. (Q we 
get: Sip + 1) °^ which is same as that 

obtained with the Jz measurement [15] . 

Let us now consider another correlated input state first 
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N N N 



FIG. 4: The variance Sip as a function of A*', the input number of photons for various states. In all the figures, the dotted 
line corresponds to the Heisenberg limit(HL) and dot-dashed to the shot-noise limit(SL). The continuous line in (a) for the 
"modified" -Yuen state, (b) for the |V')sp considered in Ref. [23| and (c) for the optimal state(for comparison S(p as obtained 
with[22] optimal POVM is also shown) 



proposed by Yuen [Tf] : 



J 1 



(18) 



Now by using Eq. (|4|, we get: (P)yu = 0. Thus the 
parity detection does not give any phase information for 
the Yuen state. The main reason for the vanishing (P)yu 
is the relative phase of 7r/2 among the two possible inputs 
at the MZI. This motivates us to consider the state with 
zero relative phase, which would be a slightly modified 
form of |V'')yu- Let us define the "modified" Yuen state 
as, 



m.yu 



E 



('^^4/2 



V2 



-1/2) I . , 



which then following Eq. (10 1 leads to: 



(19) 



(20) 



Thus using Eq. (|4| one can calculate the variance. In the 
limit of (/J — > this leads to a sub-shot noise variance as 
is shown in Fig. Qa). 



Berry and Wiseman have the so-called proposed opti- 
mal states [22], for the use of the optimal POVM pro- 
posed by Sanders and Milburn ^W.. We now consider 
this state in the context of the present work. Formally, 
the optimal state is [22] : 



(21) 



where the amplitude is given by 



1 



'~ VJ+i 



2j + 2 



It is worth noting that the optimal state is a state 
within the interferometer just like the case of the NOON 
state. So we use the operator Q as the detection operator. 
It would be a simple BS transformation to obtain the 
actual MZI input state and this result is given in Ref. [22] . 
First, we need to transform this state through the phase 



shifter 
Eq 



IV'>^|opt = 

([5| we obtain. 



e'^^^G^ I J, fi) . Then, using 



(Q) 



C^|opt 



(22) 



Plugging the above equation into Eq. ([4| we can cal- 
culate Sip. However, we have not found a closed form 
expression for Sip. Instead, Fig. Qb) shows a numeri- 
cal plot. Clearly the optimal state gives a sub-shot noise 
level and does better for large photon number. It is worth 
noting that, although one can generate the optimal state, 
the implementation of the optimal POVM [55] for which 
IV')opt designed for, requires a real time feedback. The 
variance in this case is (22]: 



Sip, 



opt 



tan 



N + 2 



N' 



(23) 



The parity detection, on the contrary, is relatively 
straightforward. The variance thus obtained due to 
optimal POVM and parity detection for the |'!/')opt 
shown in Fig. |4jb). 

Finally we consider the input state, recently suggested 
by Pezze and Smerzi 23] for achieving the HL, given by 



V2 



(24) 



The strategy employed by these authors is direct de- 
tection of number of photons at the output modes of 
MZI and applying Bayesian analysis for multiple detec- 
tions with greater confidence. Also it is important to 
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TABLE I: Phase estimate for various input states under parity detection. Here, Row 1: Coherent state at one input (fiin). 
Row 2: Number state at one input (ai„). Row 3: Dual Fock state. Row 4: Yurke state. Row 5: "Modified" Yuen state. Row 
6: State suggested by Pezze and Smerzi. Row 7: Berry- Wiseman Optimal state. Row 8: NOON state. The states in 7 and 8 
are states representing the modes a and o within the interferometer. In the Schwinger notation, Jy eigen states represent the 
internal modes a! and 6' (see Ref. 21 ). 



Input State Phase Estimate 





Fock- state Notation 


Schwinger Notation 


50 


1. 


Na|0>, 






2. 


\N)a\^\ 






3 




3 




4. 


71 [T)a T)6+ T + -'-)a T^-'-)b] 




, 1 ^72 

^ Vf (f +1) ^ 


5. 


1 r|Af+l\ |Af-l\ 1 |iV-l\ |Af+l\] 
v/2 LI 2 / a\ 2 /b^ \ 2 / a\ 2 / b\ 




sub-shot noise [Fig. 


4 


a)] 


6. 


^[|f + l)Jf-l).+ lf-l)Jf + l)J 


E. ^ l^'A^) 


sub-shot noise [Fig. 


4 


:b)] 


7. 


/ 2 Y^A^ ■ (W+2+fe)7r |Ar+fe\ |W+fe\ 
Y Af+1 Z^fe=0 ^^^^ [ 2{N+2) \\ 2 /a'\ 2 / b' 


^ .,i„(!i£+4«ll) 


sub-shot noise [Fig. 


4 


c)] 


8. 


^,\N)A^)b' + \^)a'\N)b' 




S^=j^ (HL) 



note that Eq. Q was not used to calculate the variance 
but a single interferometric measurement was used. In 
this sense, it is claimed [53] that the above state is the 
most optimal one for the HL. Such a quantum state of 
Eq. ( 24 ) was also considered in Ref |47j by Dowling with 
the measurement. 

Here we apply parity detection to the input state \ip)„^- 



Using Eq. (10), we immediately get the following result: 



(F)p. = {-iy+\di,{2^) + dUAm. (25) 

and the phase variance can be calculated using Eq. ([4]). 
We plot the result numerically Fig. |4^b), in the limit of 
(/?—!■ 0. And we clearly see the sub-shot noise limit. 
Indeed for the large number of input photons, the phase 
estimate approaches the HL. 



IV. CONCLUSIONS 

To summarize, in this article, we demonstrate the im- 
portance of the parity detection scheme in the optical 
phase estimation. By considering a combination of the 
NOON state and the dual-Fock state at the input, we 
have shown that the parity detection still gives sub-shot 



noise variance, and it reaches close to the HL for large 
input number of photons. We have also considered the 
Yuen state, Pezze and Smerzi optimal state, the Berry- 
Wiseman optimal state and have shown that we can 
achieve even smaller phase variance using the parity de- 
tection. Our results indicate that the parity detection 
acts as a unified detection scheme for precision phase 
measurements. We have summarized the results in the 
Table H 




FIG. 5: Schematic of a containing the detecting output mode 
hout from the MZI in Fig. [l] in a cavity towards measuring 
the parity of the mode as suggested in Ref. [59]. The output 
mode is reflected at mirror, M, and then fed into the cavity. 
By properly picking the curvature of the cavity morrors, one 
can trap the mode for sufficiently long time. 
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From a mathematical point of view, the parity detec- 
tion appears to be simple in comparison with any other 
detection strategy, but to experimentally realize such a 
detection scheme is not trivial. There are basically two 
different approaches to accomplish the task in laboratory. 
The first one is by employing number-resolving photodc- 
tectors gSl m 1^ [511 HI [Ml [Ml at the output de- 
tecting mode (mode bout in Fig. [l]). It should be noted 
that we do not necessarily need a photodetector at single 
photon resolution, rather a detector that would discrimi- 
nate even and odd number of photons would suffice |56j . 

Alternatively, one can measure the Wigner function. 
In Ref. [5S] it has been shown that expectation value of 
the parity is h/2 times the Wigner function. In Ref. [55] . 
it has been shown that if we can store the single mode 
field in a cavity, then we can perform the parity detection 



by measuring the Wigner function. This would require 
containing the output detecting mode bout in a cavity 
as shown in Fig. [5] By properly picking the curvature of 
the cavity morrors, one can trap the mode for sufficiently 
long time. This approach is discussed further in Ref. [60]. 
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